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We study the statistics of wave functions in a ballistic chaotic system. 

The statistical ensemble is generated by adding weak smooth random poten- 

'^ ■ tial, which allows us to apply the ballistic a-model approach. We analyze 

Q . conditions of applicability of the a-model, emphasizing the role played by the 

O I single-particle mean free path and the Lyapunov exponent due to the random 

potential. In particular, we present a resolution of the puzzle of repetitions of 
periodic orbits counted differently by the a-model and by the trace formula. 
^ ; PACS numbers: 73.21.-b, 05.45. Mt, 03.65.Sq 

in 
in 

l> ■ 1. INTRODUCTION 

o. 

The central problem in the field of quantum chaos is understanding of 
statistics of eigenfunctions and energy levels of a quantum system whose 
C^ I classical counterpart is chaotic and their relation to the underlying classical 

dynamics. The Bohigas-Giannoni-SchmilJll conjecture states that, generi- 
i-A . cally, statistical properties of levels in such a system are described (in the 

^J I leading approximation) by universal results of the random matrix theory 

(RMT). A related hypothesis concerning statistics of wave functions has 
been put forward by Berrya (see alsoEi) who conjectured that an eigenfunc- 
tion of a chaotic billiard can be represented as a random superposition of 
^ ' plane waves with a fixed absolute value k of the wave vector (determined by 



j^ ■ the energy k"^ /2m = E, where m is the mass and we set ^ = 1). This implies 



Gaussian statistics of the eigenfunction amplitude V'(r), 

V{ij} oc exp -^ f d'^rd'^r'ij*{r)C-Hr,r')^{r') 
'dedicated to Peter Wolfle on the occasion of his 60th birthday 
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determined solely by the correlation function (we consider a two-dimensional 
system) C(ri,r2) = {i/j* iri)tp{r2)) = Jo(A;|ri - raO/V". Here /3 = 1 (or 2) for 
a system with preserved (respectively broken) time reversal symmetry, V is 
the system area, and Jq the Bessel function. Note that, when taken literally, 
Eq. (H) contradicts the wave function normalization, 

d'r [(|^2(r)v,2(r')|) - {\^\r)\){\i;\r')\)] = 0, (2) 

since the integrand is equal to C^(r, r') > according to (||). Therefore, the 
limits of validity of this conjecture have to be understood. 

Despite much effort spent in this direction, no proof of these conjectures 
has been obtained via semiclassical methods, and deviations from universal- 
ity have not been calculated in a controlled way. This is because the standard 
semiclassical tool - representation of a Green's function in terms of a sum 
over classical trajectories - is only justified for times much shorter than the 
Heisenberg time tn = 27r/i/A (where A is the mean level spacing), which is 
not sufficient for the problems considered. 

On the other hand, a considerable progress has been achieved in inves- 
tigation of the statistical properties of energy levels and wave functions of 
diffusive disordered systems. In this case the supersymmetry methodD serves 
as a tool for a systematic analytical description of the level and eigenfunction 
statistics. After averaging over an ensemble of realizations of the random 
potential the problem is mapped onto a supermatrix c-model, which is fur- 
ther studied by various analytical means. This has allowed one not only to 
prove the relevance of the RMT results, but also to calculate system-specific 
deviations from universality determined by the diffusive classical dynamics, 
seeO for a review. 

This success of the diffusive cr-model gave rise to an expectation that 
the supersymmetry method may be also useful in the context of quantum 
chaos. In a seminal paper,H Muzykantskii and Khmelnitskii conjectured 
that a chaotic ballistic system can be described by a ballistic u-model, with 
the Liouville operator replacing the diffusion operator in the action. They 
presented a derivation of this model using the averaging over a white-noise 
disorder and conjectured that it remains valid in the limit of vanishing disor- 
der. Subsequently, Andreev, Agam, Simons and AltshulerQ proposed another 
derivation of the model, by considering a clean system and employing the 
energy averaging only. This led them to the conclusion that the statistical 
properties of a chaotic system can be obtained from the results found for dif- 
fusive systemsQ by replacing eigenvalues and eigenfunctions of the diffusion 
operators by those of the (properly regularized) Liouville operator. 

However, soon after publication 00 Q these conclusions were criticized 
from several points of view. PrangeH showed that the energy averaging is 
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insufficient to detect non-universal corrections to the level statistics, in view 
of the statistical noise. This indicates that one cannot study an individual 
system but should rather consider some ensemble of systems. The necessity 
of an ensemble averaging has been also pointed out by Zirnbauercl who con- 
centrated on rigorous formulation and justification of the Bohigas-Giannoni- 
Schmit conjecture of universality. Bogomolny and Keatingly emphasized the 
discrepancy between the prediction of the ballistic cr-model for the smooth 
part of the level correlation function and the result of the diagonal approxi- 
mation to the semiclassical trace formula. Specifically, repetitions of periodic 
orbits are counted differently within the two approaches. These two issues 
(ensemble averaging and repetitions) point to serious problems with the bal- 
listic (T-model approach and require clarification. This will be done below. 



2. BALLISTIC cr-MODEL FROM AVERAGING OVER 
SMOOTH DISORDER 

Let us come back to the derivation of the cr-model. It turns out that 
the derivations proposed in both Refso andQ are, in fact, not quite correct. 
Specifically, if averaging over a white-noise disorder with a mean free path 
I is used (as inD), the ballistic a-model is only valid for momenta q <^ l~^ , 
i.e. in the region of validity of the diffusive cr-model, but not in the ballistic 
range of larger momenta.O On the other hand, the energy averaging of Ref.l3 
leaves a continuum of zero modes allowing for arbitrary fluctuations of the 
supermatrix field transverse to the energy shell and spoiling the derivation 
of the ballistic cr-model, which reouires that these fiuctuations be frozen.Q 

We note that the jiiethods oia anclEI may be considered as opposite ex- 
tremes: while in Ref.Q the averaging over a random potential with zero 
correlation length d = was proposed, the energy averaging of Ref.Q cor- 
responds to a random potential with d = oo. As often happens, the truth 
lies in between: the ballistic cr-model can be obtained if one averages over a 
smooth random potential with a finite correlation length d. In fact, thistype 
of derivation had been performed for the first time by Wolfle and Bhattll3 ten 
years before the notion of the ballistic a-model was introduced. However, 
the aim of Ref.t3 was demonstration of the applicability of the diffusive cr- 
model to the problem of a smooth random potential, and the ballistic action 
appeared only implicitly as an intermediate step of the calculation. This 
derivation was generalized to the case of a random magnetic field inJ13 yet 
again the main point of that work was the diffusive action on distances ex- 
ceeding the transport mean free path. More recently, it was emphasizedHaO 
that averaging over a smooth disorder is exactly the proper way of derivation 
of the ballistic cr-model. This established the connection between Refs.113113 
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on the one hand, and Refs-BQ on the other hand. 

Starting from a system with a Hamiltonian H, we generate a statistical 
ensembleEJ by adding a random potential U{v) characterized by a correlation 
function W{r — r') = {U{r)U{r')) with a correlation length d. Parameters of 
this random potential are assumed to satisfy k^^ '^ d ^ Ig ^ Itj-, where Is = 
VTs, kr = VTtj., V is the velocity, and rs(rtr) is the single-particle (respectively 
transport) relaxation time. Note that the potential is smooth, kd ^ 1, so 
that kr/ls ~ {kd)"^ » 1. In order to apply this idea to a ballistic system 
of a characteristic size L, we require Ig <^ L <^ Itr- The condition Ztr ^ L 
preserves the ballistic nature of the system, while the inequality Is -^ L 
guarantees that the ensemble of quantum systems is large enough to produce 
meaningful result. 

After the ensemble averaging, the problem can be reducedEJ't3't3'll3 to 
a "non-local ballistic cr-model" of a supermatrix field Q(r, n) with the action 
(for definiteness, we consider the case /3 = 2) 



S[Q] = Strln 



■KV 
T 



E-H + ujA-- fdn'Qir, n)w{n, n'] 
I dVdndn'Str Q{t, n)u;(n, n')Q(r, n'), (3) 

V is the density of states, w(n, n') = InuW^kln — n'\) is the scattering 
cross-section for the random potential, and n is a unit vector characterizing 
the direction of velocity on the energy surface. Note that despite the non- 
local Str In form, the action (^ is not an exact representation of the original 
problem, but rather a low-energy theory with only soft modes kept: the 
momentum variable of Q is constrained to the energy surface, and Q satisfies 
the usual u-model condition Q^ = 1. To obtain the ballistic cr-model in the 
local form (as proposed irJj'Q), one has to perform a gradient- and frequency- 
expansion, which is justified provided qlg ^ 1, coTg <C 1. The result iglflO 

S[Q] = TTiy f d^rdnStr AT-\r,n)CT{r,n) - ^-^AQir,n) 

+ ^ / d^rdndn'StrQ(r,n)u;(n,n')(5(r,n'), (4) 

where Q{r,n) = T{r,n)AT~^{r,n). The symbol C denotes the Liouville 
operator characterizing the classical motion; for a billiard C = vnV supple- 
mented by appropriate boundary conditions .113 The local ballistic cr-model 
(Q) is thus only applicable on length scales ^ Is- From the point of view of 
the semiclassical (periodic-orbit) theory, this corresponds to the condition 
of applicability of the diagonal approximation. On shorter scales, one has to 
use the more general, non-local, form (^). 
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3. EIGENFUNCTION STATISTICS IN A BALLISTIC 

SYSTEM 

The two-point correlation function of the wave function intensities is 
expressed in this approach asErlij 

(|V2(r)V'2(r')|) = hm^([Gii(r,r)G22(r',r')+Gi2(r,r')G2i(r',r)])5[Q], 
77-^0 vr 

where G is the Green's function in the field Q, 



G 



E - H + irjK - - dn'Q{r, n')w{n, n') 



-1 



(6) 



and the subscripts 1,2 refer to the advanced-retarded decomposition (the 
boson-boson components being implied). Here t] is an infinitesimal positive, 
and CO in the action S[Q] is given by a; = irj. 

We first evaluate Eq. (||) in the zero-mode approximation, Q{r) = Qq. 
The Green's function (y) is given in the leading order by 

Go(r,r') = iImGR{r,r')Qo + ReGR{r,r'), (7) 

GRir,r') = {r\iE-H + i/2rs)-'\r'). (8) 

Substituting Eq. (0) in (g) and expanding the action (y) up to the linear-in- 
?7 term, S[Q] ~ vrz^j^yStrQoA, one finds, in a full analogy with the case of 
diffusive systems, Ej 

V^{\iJ'{rl)^//{r2)\)c^l + kg{rl,r2)■, (9) 

kg{r, r') = ImGij(r, r')ImGfl(r', r)/{7riyf, (10) 

with the two contributions on the r.h.s. of (^ originating from the terms 
(G11G22) and (G12G21) in @, respectively. The result @, corresponding 
exactly to the conjecture (|l|) of the Gaussian statistics, is in conflict with 
the wave function normalization, as explained above. 

To resolve this problem, we evaluate the term (G11G22) more accurately 
by expanding the Green's function (^) to the order ry and the action (^ to the 
order rj'^. While these terms (usually neglected in the c- model calculations) 
are of the next order in ry and may be naively thought to vanish in the limit 
r) —)■ performed in @, this is not so, since Qo oc r/~-^. As a result, we get 
in the zero-mode approximation 

VWij\ri)ip\r2)\)zM - 1 = kg{ri,r2) - k,{vi) - ~kg{r2) + kg (11) 

(terms of still higher orders in 77 produce corrections small in the parameter 
At, < 1), where kq{r) = V-^ f d^r'kg{r, r'), kq = V-^ J £rd^r'kg{r, r'). 
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The contribution of non-zero modes is expressed in terms of the cr-model 
propagator ,0 

V^{\iI,^{vi)^^{v2)\)nzm = nB(ri, rs), (12) 

where riB(ri,r2) = nB(ri,r2) — Ilg (ri,r2) describes the (integrated over 
direction of velocity) probabihty of classical propagation from ri to r2, 

nB(ri,r2) = / / dnidn2P(rini,r2n2), 

CV = (vri/)-! [6{yi - r2)5(ni - n2) - V'^] , (13) 

with the contribution 11^ (ri, r2)pQf direct propagation (before the first event 
of disorder scattering) excluded JIZI 

We analyze now the total result given by the sum of (O) and ([l2|). First 
of all, we stress that it satisfies exactly the condition (||) of wave function 
normalization. Next, we consider sufficiently short distances, |ri — r2| <^ h- 
In this case the correlation function is dominated by the first term in the 
r.h.s. of Eq. (0), returning us to the result (^). Furthermore, we can 
generalize this result to higher correlation functions, 

{r{ri)ijir[)...rirnMrn)) 



2V(n-l)l 



X lim(2vr.,)"-i (E IT ^^-P.cJr^'M.)) ) > (14) 

where the summation goes over all permutations a of the set {1, 2, . . . , n}. 
Pi = \ for i = 1, . . . , n — 1, and pn = 2. If all the points r^, r[ are within a 
distance <C Is from each other, the leading contribution to this correlation 
function is given by the zero-mode approximation with higher-order terms in 
r] neglected [i.e. by the same approximation which leads to Eq. @j\, yielding 



y"(r(ri)V(r'i)...V'*(r„)V^(r:,)) = 5;n/^(r-<w)' 
/i.(r,r') = -ImGfi(r,r')/(^z.). (15) 



This result is identical to the statement of the Gaussian statistics of eigen- 
functions conjectured inoQ We have thus proven that Eq. (|l|) holds within a 
spatial region of an extension <S Ig-, with the kernel C(ri, r2) = /^(ri, r2)/y 
given by Eq. (|lq). 

We turn now to the behavior of the correlator (|V^^(ri)'(/'^(r2)|) at larger 
separations |ri— r2l ^ Is- In this situation, the correlations are dominated by 
the contribution ([l^ ) of non-zero modes. Let us further note that the smooth 
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part of the zero- mode contribution (11) (i.e. with Friedel-type oscillations 



neglected) is exactly equal to H^ . Therefore, the smoothed correlation 
function is given by the classical propagator, 

VWi;\ri)i:\r2)\)smootl. - 1 = nB(ri,r2), (16) 

independent of the relation between |ri — r2| and Ig. The mean free path 
Is manifests itself only in setting the scale on which the oscillatory part of 
(|-i/'^(ri)V'^(r2)|) gets damped. 

To summarize, the disorder averaging generates the scale Is which sep- 
arates the regions of applicability of the Gaussian statistics and of the qua- 
siclassical theory (Q). However, this is not the full story yet. As we demon- 
strate below, the disorder averaging induces one more scale which plays a 
crucial role for the problem considered. 

4. PROBLEM OF REPETITIONS 

As was mentioned in Sec. |I], Ref.ll3 emphasized discrepancy between the 
cr-model and the trace formula as concerned the counting of repetitions in 
the expression for the level correlation function. As we demonstrate in this 
Section, this problem, while not affecting (except for a non-generic situation 
when the points ri and r2 lie close to one and the same short periodic orbit) 
the lowest-order correlation function (|^^(ri)^^(r2)p) considered in Sec. y, 
is also relevant to higher-order correlators of wave function amplitudes. We 
will clarify the origin of this puzzling discrepancy and formulate the condi- 
tions under which each of the two results apply. We assume /? = 1 in this 
section (so that eigenfunctions are real) and start our consideration from the 
correlation function 

rW(r - r') = {ij^{r)Mr)M^)Mr)M^')Mr')M^')Mr')), (17) 

where ip^^u,p,a are different eigenfunctions with sufficiently close energies. 
We will further assume that Is <^ |r — r'| <^ Itj-, i.e. the distance |r — r'| 
is in the expected range of applicability of the ballistic cr-model (^). The 
objects of the type (|l7|) naturally arise when one studies fluctuations of 
matrix elements of the electron-electron interaction, which are important 
for statistics of electron transport through quantum dots, see.ll3Ej 

We first consider the situation of a sufficiently large system size L {e.g., 
we can assume a diffusive system, L ^> Itr); in this case the value of L will 
be irrelevant for the results (except for normalization of wave functions). 
In the end of the section we will generalize our conclusions to the case of 
a ballistic system. We start by performing an averaging over an auxiliary 
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random potential with a single-particle mean free path Ig = vfg satisfying 
\r — r'\ <^ Is <^ L (which is much weaker than our "main" random potential 
and will thus not enter the final result). This auxiliary averaging allows us 
to present (^) in the form 

r(4)(r-r') = (v^^.)-4([ImGK(r,r')]^) = ^(7ri.)-^G|(r, r')Gi(r,r')), (18) 

where Gji^a = {E — H ± i/2fs)~^ is a Green's function in a given realiza- 
tion of the "main" random potential, the averaging over which remains to 
be performed in the r.h.s. of (18), as indicated by (. . .). The cumulant 



((V'^(r)'0^(i''))) of ^ single eigenfunction can be also reduced to the form 
(p^) , if one uses the cr-model approach to average over the auxiliary random 
potential, in full analogy with the derivation of Eqs. (P), ( [lO| ) (in view of 
|r — r'l <^ Ig the zero-mode approximaticp^ is appropriate). An analogous 
trick of an auxiliary averaging was used inc3 in order to study ballistic wave 
function correlations in a random magnetic field. 

Let us first evaluate the r.h.s. of ( p!8| ) using the ballistic a-model. The 
calculation is trivial and yields 

(2vr2r.2)-2(G|(r,r')Gi(r,r')) = 2n|(r, r') =. ^^^^^^_ ^,^^, ■ (19) 



The result (|T9|) is fully transparent from the point of view of diagrammat- 
ics: there are two possibilities to couple the retarded and advanced Green's 
functions in two "ballistic-diffuson" ladders, yielding the factor of 2 in front 
of the squared propagator n^(r,r'). 

We are going to show now that the result (^) of the ballistic cr-model 
is only correct for sufficiently large distances, |r — r'[ ^ II (the scale II 
will be specified below), while in the opposite limit, |r — r'| <C II, it is 
wrong byiactor of 2. To evaluate the r.h.s. of (^), we use the path integral 
approach£ll The product of the four Green's functions in ([iq ) can be written 
as (we set r = and R = r' — r) 



4 



foo pri{ti)='R 



roo i-ri(ti)=in. 

{Gl{0,K)G\{0,K)) =11 dU Vr, exp[i5kin - ^w] 

i=l-JO Jri{0)=0 

X exp[z^(ti +t2-t3- U) - {h + t2 + t3 + t4)/rs]; (20) 

5kin = y n ' dti-l + f ' dti-l - f ' dtvl - I dti-lj ; 

'S'W = ■t(5'ii + 5*22 + Ss3 + 544) -I- 512 + 5*34 — 5i3 — Si4 — S23 — 824; 

S^j= /" ' / ' W{ri{t) - rj{t'))dtdt', (21) 

Jo Jo 
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where the paths ri(ti),r2(t2) correspond to the retarded, and r3(i3),r4(t4) 
to the advanced Green's functions. It is useful to perform the change of 
variables, introducing T = (ti+i2+i3+t4)/4, t_ = (ti+t2— is — i4)/2, T12 = 
h - t2, T34 = ts - t4, and R+ = (ri + r2 + rs + r4)/4, R„ = (ri + r2 - rs - 
r4)/2, ri2 = ri — r2, r34 = r3 — r4. The kinetic part of the action then reads 

Since we are interested in the ballistic scales (<C hr) it is convenient to 
split r, ri2,r34 into components parallel (||) and perpendicular (_L) to R. As 
shown below, S'w depends only on the _L-components. Shifting the parallel 
components of ri2,r34 via pi2 = r^g — -R+T12/T', P34 = r^^ — R^Tu/T^ we can 
perform integration over P12, P3'i,Ti2, and T34, which yields the factor (T/R)'^. 
Furthermore, integration over i?i[ , i?_ , t_ , and T produces the factor m/2E 
and sets T = R/v with v = (2E jm)^''^ . Finally, introducing 

p± = (rf2 ± r3^4)/2, r± = i?i, (23) 

we simplify the action to the form 

^kin = m dt {2r+r_ + p+p_} . (24) 

Jo 

Now we turn to the disorder-induced part of the action, S'w (where all 
ti can be approximated by T, seeE^). Expanding the correlator W up to the 
second order in r, we obtain 

T rT 

2 



S^^ U{p-{t),p+{t))dt-2 {G{p^{t)) + G{p+{t))]rt{t)dt, (25) 

JO JQ 

where 

U{y, y') = 2{F{y) + F{y')] - F{y + y') - F{y - y'). (26) 

Here the functions F[y) and G{y) are given by 

/■oo J ^00 J q2 

F{y)^ —[W{x,Q)-W{x,y)l G{y) ^ — w{x,y), (27) 

JQ V Jo V dy^ 

with the following asymptotic behavior F{y ^ d) ~ — G(0)y^/2, F{y ^ 
d) ~ T.-i, G(0) = -m^wVrtr, and G(y » d) ^ 0. 

We examine first the contribution of the region p+ , P- ^ d, where all 
four paths are uncorrelated. Then the action (|25|) is large, 5w — 2TF(y ^ 
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d) = 2T/ts, since the phases acquired by the waves travehng along these 
paths do not cancel each other. Thus the contribution of this region decays 
exponentially at R ^ Ig (which is the scale of our interest), in the same 
way as the single Green's function does. We conclude that the path integral 
(|20| ) is dominated by paths with at least one of p± much less than d, which 
corresponds to correlations between "retarded" and "advanced" paths. 

Let us now consider the contribution of /9+ ^ d, p_ <C d, when the 
two "ballistic diffusons" (formed by the pairs of paths 1,3 and 2,4) are not 
correlated. Then Sw acquires a form 



Sw 



V V 



dt, 



(28) 



where wq = —vG{0) ~ v{d'^Ts)~^. Combining Eqs. ( p^ ) and (28) we see that 
the variables r± and fl± separate, each of the two pairs being characterized 
by the diffusive actioncil (describing the dynamics of the "center of mass" r+ 
and the distance between the two "diffusons" /9+, respectively). This gives 
squared propagator n^(r,r'), as in Eq. (p^). The region p_ S> d, p+ <^ d 
(diffusons formed by paths 1,4 and 2,3) yields an identical contribution, 
reproducing the factor 2 in the a- model result ([T9|). 

However, the approximation ( p8| ) is only justified at sufficiently large T. 
Indeed, the action (|^ 



(28) implies a diffusion of the velocity p+fv with 
the diffusion coefficient l/rtr. Taking into account the boundary conditions 
/3_l_(0) = p+{T) = 0, it requires a time of the order of 



TL ~ Ttrid/l 



tr 



l2/3 



(29) 



for p+ to reach a value ~ d, which implies T ^ r^ as a condition of validity 
of Eq. (28) and thus of the result (19). As we discuss below, the scale tl 
has a meaning of the inverse Lyapunov exponent; hence the subscript "L". 
At shorter distances R <^Il = vtl , both pj^ and p- are small compared 
to d. Expanding (p5[) up to the second order in both /9+ and p- we have 



S'w 



—r_{t) + —p^{t)p_{t) 



dt, 



(30) 



where W2 = —vG"{0)/2 ~ wo/d^ and we neglected cross-terms of the type 



rj^/9^, which are small compared to the first term in the integrand in (pO|). We 
note that the variables r± and p± again separate. The first pair corresponds 
to the diffusion of r+ [with a diffusion coefficient two times smaller than in 
Eq. ([2^)], while the second one now describes the divergence of the two pairs 
of paths. The differential equation for a Green's function which corresponds 
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to the /9-part of the action reads 

and we are interested in the quantity g{0,0,t), the correlation function (19) 



(2vr2z.2)-2(G|(r,r')Gi(r,r')) = ^^^g{0,0,R/v) . (32) 



being given by 

2 

Now we analyze the solution of (|3l|). Clearly, tl = (m^t;/ti;2)^''^ sets 
a characteristic time scale for this equation. It is remarkable that upon 
the Fourier transformation, />+ —^ i{mv)^^d/dcl) and d/dp+ —f —imvcj), the 
operator in thel.h.s. of (31) coincides with the one obtained by Aleiner and 



Larkin in Ref.l^ from the analysis of divergence of classical paths in a weak 
smooth random potential. As shown in,E3 the scale tl has the meaning of 
a corresponding inverse Lyapunov exponent. We have thus demonstrated 
that the same scale arises in the treatment of a quantum (Is ^ d) random 
potential, despite the diffractive nature of the scattering.lH3 

For t <^tl one can neglect the term p\p'^_ in the l.h.s of (|3l|), yielding 



m ip+P~ 

5((/)+,p_,t) ~ -— -exp — , t<.TL. (33) 

ZTit mt 

In fact, the solution remains a function of the product p+p- for all t, and 
Eq. (^) can be reduced, in dimensionless variables r = —i^'^t/rL, r = 
(z/4)^^'^(mp_|_p_/T/^)^'^, to an equation for a 2D anharmonic oscillator, 

{-idr - r-^drvdr + r^)g{r, r) = 6^{r)6{T). (34) 

Therefore, at r ^ 1 (corresponding to t ^ tl) its solution decays expo- 
nentially, g oc exp{— z^'^eot/Ti}, where eg is the (dimensionless) oscillator 
ground-state energy. However, as explained above, at this time scales the 
approximation p+,p- <C d leading to Eq. (^) loses its validity, and the 
result is dominated by the contributions with either p^ ^ d or p^ ^ d 
described by the action ( p8| ) and yielding the result (|T9|). 
We thus conclude that 

(2^2^2)-2(G|(r,r')Gi(r,rO) = jJ]^HR/h), (35) 

where J^{x) is a parameterless function with asymptotics T{x ^ 1) = 1 and 
J^{x » 1) = 2. To calculate the crossover function T{x), one has to solve a 
differential equationE^ 

I " i^dS^^^^P^'^-^) 9(P+^P-^t) = mS{p+)Sip-), (36) 
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with hl{p^, P-) defined in (pq); the desired correlation function is then given 
by Eq. (pi]). The sokition in the crossover region depends on the specific 
form of the random potential correlation function W{r) and, in general, can 
only be found numerically; however, the explicit form of the function J-{x) 
is not important for our analysis. A similar analysis can be performed for 
higher-order correlation function (G^(r, r')G^(r,r')), yielding the crossover 
from ^(„)(x <C 1) = 1 to ^(n)(x ^ 1) = n\. 

Therefore, the ballistic o"- model result is only valid at distances R^ II- 
At shorter scales, the four paths do not split into two pairs yielding each a 
ballistic diffuson (as the cr-model approach assumes) but rather propagate 
together, remaining all strongly correlated. In the diagrammatic language 
this means that all the four Green's functions are coupled by impurity lines. 
This is closely related to the non-Markovian memory effectsE3 showing up 
in transport in a smooth random potential. 

For ballistic systems this implies that the cr-model approach is only 
valid provided the random potential over which the averaging is performed 
is sufficiently strong, so that II <^ L. Note that such an averaging can 
be characterized as "strongly invasive": the Lyapunov exponent due to the 
random potential t^ is much larger than that of the clean system itself 
(which is ~ v/L for a generic system of a size L). Therefore, although such 
a random potential does not essentially influence (in view of /tr ^ L) the 
positions of Ruelle resonances governing the relaxation rate in the chaotic 
system, it strongly affects the Lyapunov exponents. 

In the opposite case II ^ L (corresponding to a "non-invasive" aver- 
aging, keeping intact all the classical features of the system) the ballistic cr- 
model does not give correct results for non-universal level and eigenfunction 
correlations. More precisely, it is still valid when applied to the lowest-order 



correlation function {\^ {r)'ip {r')\), see Eq. (16). However, when one con- 
siders higher-order correlation functions (or, in fact, (|^^(r)^'^(r')|) with r 
and r' located close to a short periodic orbit), the approach breaks down, 
since it does not take into account correlations between four or more close 
paths. This is precisely the problem which shows up in the counting of repe- 
titions in the formula for the level correlation function. In that case, relevant 
paths for both Gr and Ga wind n times in the vicinity of a periodic orbit. 
If L <C II, all these paths are correlated (in the same way as four Green's 
functions in the example considered above). The o"- model neglects this fact, 
combining the two paths ("R" and "A") in a ladder, which can be done in n 
possible ways. This leads to an overestimate of the contribution by a factor 
of n (similarly to the overestimate of the correlation function (|l9| ) at i? <C /l 
by a factor of 2). Therefore, the results of the cr-modelQ and of the trace 
formulcJiy for the smooth part of the level correlation function correspond to 
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different types of averaging: the former to II <C L, while the latter to II ^ L 
(in both the cases the condition Is <^ L should be assumed, see Sec. ^). 

It is worth noting that we assumed the clean system to be characterized 
by a single scale L, so that the corresponding Lyapunov exponent is ~ L~^. 
An example of a system which does not belong to this class is given by 
billiard with a surface roughness producing a diffuse boundary scattering.! 
In this case repetitions eive a negligibly small contribution, and the a-model 
approach as applied intll is justified with the only assumption Is <^ L. 

Finally, let us note a close similarity between the condition II > L oi 
the failure of the a-model and the condition //, > Re (with Re the cyclotron 
radius) under which memory effects develop into adiabaticity which affects 
strongly magnetoresistivity in a smooth disorderH3 In both cases, correla- 
tions between multiple traversals of a periodic (respectively cyclotron) orbit 
lead to breakdown of the Boltzmann kinetic equation. 



5. CONCLUSIONS 

We have studied the statistics of wave function in a chaotic system on 
the ballistic scale. To define the statistical ensemble, we have added to 
the system a smooth quantum random potential with a correlation length 
d ^ k~^. We required the corresponding quantum (single-particle) mean 
free path to satisfy Ig ^ L, which ensures that the ensemble of quantum sys- 
tems is sufficiently large and provides meaningful statistics. Using the ballis- 
tic (T-model approach, we have shown that on scales <C h the waye function 
statistics is Gaussian, Eq. (Il5|), proving the Berry's conjecture. Bh On larger 
scales Friedel-type oscillations in the correlation function {\ip'^ [r)ip'^ [r')\) are 
exponentially damped, and the latter is given by the classical ballistic prop- 



agator, Eq. (12). However, an attempt to use the u-model approach for 
calculation of higher-order correlation functions in this regime makes us to 
face a problem: the cr-model does not take into account strong correlations 
between 4 (or more) Green's functions dominated by paths which are all 
close to each other. Considering the correlation function ([iq ) as an exam- 
ple and performing its detailed analysis, we demonstrated that the ballistic 
(J- model result is only correct on distances R ^ II = vtl, where r^ is 
the classical Lyapunov exponent associated with the random potential, see 
Eq. (P9|) . It is remarkable that the Lyapunov scale II arises despite the fact 
that the random potential is assumed to be quantum, with Ig ^ d. On 
shorter distances, R <^Il, the u-model result for (|l8|) is wrong by factor of 
two. In this regime the path integral approach is more appropriate; it allows 
one also to describe (at least, in principle) the crossover between the two 
regimes {R ~ II)- 
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A completely analogous situation is encountered when one studies the 
level correlation function. Ourresults thus resolve the discrepancy between 
the predictions of the cr-modeO and of the trace formulali3 for the smooth 
part of the two-level correlation function. The cr-model result is generically 
valid provided II *C L, where L is the system size. This condition means that 
the averaging is "invasive" , since it strongly affects the Lyapunov exponent 
(although it does not affect, in view of the assumption Itr ^ L, positions of 
the Ruelle resonances ).E3 In the opposite case of a "non-invasive" (preserv- 
ing all classical parameters of the system) averaging, l^ ^ L, the cr-model 
approach neglecting correlations between multiple traversals of a periodic 
orbit loses its validity, and the smooth part of the level correlation function 
is given by the trace formula.t3 Calculation of the system-specific oscillatory 
contribution in this situation remains an open problem. Such contributions 
are related to the behavior of the spectral form-factor at times close to the 
Heisenberg time tn, which is beyond the region of validity of the trace for- 
mula. Bogomolny and Keatingl2l proposed a procedure allowing to obtain 
the oscillatory contribution from the trace formula. However, in view of 
the ad hoc nature of their suggestion, the status of the result is unclear. 
On the other hand, the cr-model approach, which would be more appropri- 
ate for an analysis of contributions of this kind (non-perturbative from the 
field-theoretical point of view), does not treat properly correlations of mul- 
tiple close paths. A natural idea would be to combine the two approaches, 
i.e. to formulate a field theory which would contain not only conventional 
"2-diffusons" (i.e. those generated by a product GrGa of two Green's func- 
tions) P(rn, r'n'), but also "4-diffusons" , "6-diffusons" , etc. This remains a 
challenge for future research. 
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